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Abstract. Electrostatic oscillations in cold electron- 
positron plasmas can be coupled to a propagating elec- 
tromagnetic mode if the background magnetic field is in- 
homogencous. Previous work considered this coupling in 
the quasi-linear regime, successfully simulating the elec- 
tromagnetic mode. Here we present a stability analysis of 
the non-linear problem, perturbed from dynamical equi- 
librium, in order to gain some insight into the modes 
present in the system. 



1. Introduction 

Pulsar magnetospheres are composed of magnetised 
electron-positron plasmas. Conventional modelling uses 
the fact that the dominant Lorentz force produces a host 
of relativistic charged particles, each of which radiates 
strongly and stochastically, producing 7-ray photons. Such 
single-particle models have been explored as possible radi- 
ation sources (da Costa & Kahn [TWf| . but have not been 
able to recover the highest energy radiation observed. 

However this strict single particle approach may be 
challenged as there are physical processes which develop 
in the rest frame of the plasma, and are strongly ampli- 
fied when translated in the laboratory frame of the pulsar. 
In cold non-relativistic plasma theory non-linear electro- 
static oscillations in electron-positron plasmas develop a 
density instability in which the density of both species 
grows sharply at the edges of the oscillation site (da Costa 
et al. I200ip . Folding thermodynamics into the system pro- 
vides a possible mechanism for avoiding the onset of the 
instability, since pressure effects would oppose the den- 
sity build up. However, a more advantageous mechanism 
would be coupling the oscillation to an electromagnetic 
mode via an inhomogeneous background magnetic field. 
This would allow energy to be radiated away from the 
oscillation, quenching the density instability and giving 
a source of radiation in the pulsar magnctosphere. Pre- 
vious work has studied this mechanism and successfully 
simulated the coupling between the oscillation and the 
electromagnetic mode in the quasi-linear regime (Diver 
et al. I2002|) . Therefore the study of these processes is of 



the greatest importance for the dynamics of the magneto- 
sphere, and the radiation mechanisms of pulsars. 

This paper addresses the basic mathematical formula- 
tion of the electron-positron cold magnctoplasma and the 
mode coupling mechanism in Section 2. In Section 3 lin- 
ear analysis of the equations perturbed from dynamical 
equilibrium is performed to gain an understanding of the 
modes present in the plasma and their stability. 



2. Model equations 

In cylindrical polar coordinates (r, 9, z) with an axial mag- 
netic field B — zB z and an electric field in the r, 6 plane 
~E r fi the full non-linear model equations for a cold electron- 
positron plasma are 



rh + + (rn+u r ) = 





(1) 


rh— + (rri-V r ) = 





(2) 


u r + u r u' r — Ug/r = 


(e/m)(E r + u e B z ) 


(3) 


ue + u r u' r + u r ug/r = 


(e/m)(Eg - u r B z ) 


(4) 


v r + v r v' r — Vg/r = 


-(e/m)(E r + vgB z ) 


(5) 


V0 + v r v' r + v r vg/r = 


-(e/m)(Eg-v r B z ) 


(0) 


(rE r )' = 


(e/e )r(n + - n_) 


(7) 


= 


—E r /c 2 — noe(n + u r 


— n-V r ) (8) 


(rEg)' = 


-rB z 


(9) 


K = 


-Eg/c 2 - /j, e(n + u e 


-n.vg) (10) 



where n + , n_ are the positron and electron number densi- 
ties; u, v are the positron and electron velocities; ' denotes 
d/dt; ' denotes d/dr; ([7]) is Poisson's equation; ^ is the 
single z component of the induction equation; and (1101) . 
|8]) arc the 8 and r components of the V x B equation. 
The equations can be recast into a form that highlights 
the symmetry of the electron-positron plasma, this can be 
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done via the following relations 



n = 


— ( T\,a_ 4-77, \ £ / T),r\ 

2 v fc + * it— } t,/ 


(ID 


A 


-(n+ -n-)£/n 


(12) 


a = 


— (Ur + Vr) /cOriL 

2 v ' > / 1 u 


(13) 


S = 


i 

- (U r - t>r) /Woi 

z 


(14) 


X = 


i (u e + we) /w i 


(15) 


C = 


1 

- (ue - we) A->o^ 


(lb) 






(17) 


p = 




e = 




(18) 




ITIUJqL 


A) + P = 


eB z 
mujQ 


(19) 


r = 


£L 


(20) 


t = 


t/uj 


(21) 



where L,l/u>o are characteristic length, time scales and 
Po, (3 are the equilibrium and perturbed magnetic held 
respectively. The system of governing equations then be- 



comes 








t = 


- (Eo- + AS)' 




(22) 


A = 


-(Aa + ES)' 




(23) 


a = 


-i(a 2 +<5 2 )' + ( X 2 + C 2 )/e + C(/?o + 


P) 


(24) 


5 = 


-(a5)' + 2 X C/Z + P + xWo + P) 




(25) 


X = 


-x'a-t'S-ixa + CS) /£-5(p o +0) 




(26) 


C - 


- 6yd - K + S X ) It + - o (ft, - 




(27) 


p = 


-|(Atr + E<y) 




(28) 


6 = 


-pC9b+j9)'-|(SC + Ax) 




(29) 


$ = 


-6' - 9/Z 




(30) 



where p = c 2 /(w 2 L 2 ). 

2.1. The coupling mechanism 

Consider a homogeneous background magnetic field. 
When the density of the plasma is perturbed, resulting in a 
charge imbalance, a radial electric field is created that ac- 
celerates the electrons and positrons in radially opposite 
directions. The plasma collectively is responding to the 
presence of the electric field and is trying to restore equi- 
librium. The magnetic field causes the particle trajectories 
of both species to participate in partial Larmor orbits with 
the same azimuthal velocity. The particles overshoot their 
initial positions, due to their acquired kinetic energy, and 



produce a new charge imbalance that they again try to 
correct. This is one half period of an electrostatic oscilla- 
tion. 

Introducing an inhomogeneous magnetic field induces 
a B x VB drift which causes a net current density in 
the azimuthal motion of the plasma during the electro- 
static oscillation. The resulting current density induces ax- 
ial magnetic field fluctuations which propagate away from 
the electrostatic oscillation site. This mechanism causes 
the electrostatic oscillations to couple to a propagating 
electromagnetic mode in the plasma. 

This effect has been fully investigated in the quasi- 
linear regime, in which the variables were perturbed from 
a uniform equilibrium. However having an inhomogeneous 
magnetic field permeating the plasma requires the equi- 
librium to be non-uniform for the system to be self- 
consistent. 

2.2. Dynamical equilibrium 

In the equilibrium situation d/dr = and A = p = 9 = 
(3 = 0. The resulting dynamical equilibrium initial condi- 
tion is described by 



Co = -y-A) (31) 

= -Co 2 (33) 
S = kq/o-q (34) 



where subscripts denote equilibrium value and kq,ki 
are constants. The equilibrium equations describe the self- 
consistent response of the plasma to a prescribed back- 
ground inhomogeneous magnetic field, (3q. Conversely they 
describe how the plasma must behave in order to generate 
the same magnetic field. Equation ([33]) defines the kinetic 
energy conservation of the equilibrium flow; equation (I34|) 
defines the conservation of the total number density flux; 
and equations (|32|) .([3"T j) describe the magnetic field gener- 
ation via the (q, (3q coupling. 

Linearise the full set of governing equations and look 
at large values of £. In this regime the magnetic field 
is asymptotically tending to a constant value, (3q = 
constant, corresponding to no motion of the plasma in 
the azimuthal direction Co = Xo = 0, saying that radial 
flow dominates. This requires <To = constant to be consis- 
tent with equation ([33)) and implies that in this parameter 
set that there is a net motion of the plasma in the radial 
direction with no net current, 5q = 0. The governing equa- 
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tions then become, 



v 


— {L, <T + ZjCTq) 




(6b) 


A = 


-(<t A + 5T, )' 




(36) 


(7 = 


—una' + BnC 




(37) 


5 = 


-a a 5' + p + fox 




(38) 


X = 


-a x' - 0oX/£ - 


■Spo 


(39) 


C = 


-<T0(' - ^C/C - 


<j P - p a a + 9 


(40) 


p = 


2 

--(Aa + S (5) 




(41) 


9 = 


- P f3' - 2S C/e 




(42) 


P = 


-9' - e/i 




(43) 



3. Stability analysis 

Upon inspection of the linearised governing equations it 
is evident that they can be split into two independent, 
self-consistent sets namely equations (|36|) , (|38t . (|39)) , (|4T|) 
forming one set and equations (|55 ]) .([3T |l .(|lD ll .(|i2 ]) and ([33"]) 
forming the other. This simplification was exploited to ob- 
tain the following solutions. 

3.1. Electrostatic solution 

The electrostatic solution is characterised by (3 = Po, 9 = 
and C = 0. Substituting these conditions into the govern- 
ing equations yields, 



A = -(a A + <5So)' 
5 = -o-qS' + p + PoX 
X = -oox' - o-qx/S, - PoS 
p = -2(Aa Q + E <5)/£ 



(44) 
(45) 
(46) 
(47) 



Combing equations ([4"5")) and (T4"6")) to eliminate x and (HU) 
and (T4"4")) to eliminate A produces two differential equa- 
tions both involving S and p. Substituting one differential 
expression into the other to eliminate p yields the partial 
differential equation 



S + 2a S' 



2 c/l 

a 5 



o-oS/Z + aid' /t + {Pi + 2S /0 5 

-2 ( 7 C 1 /e = (48) 

and setting S — 



where C\ is a constant. Setting C\ 
?/(£) exp (— iojt) yields 



<V" + (*o/£ - 2iu<r )y' + {pi + 2£ /£ - uj 2 

-iua /£)y = 



(49) 



This has complete solution 



S (£, t) = f-V3 e *-«/^-r) [C 2 M(- 



+C 3 W{- 



iS Q 2iM, 
2iPoH 



o-oPo 
iS 



,0. 



)] (50) 



where M, W are Whittaker functions of the first and sec- 
ond kind respectively and C3 are constants. M and W 



are related to the confluent hypergeometric functions M. 
and U as follows 

M(k,(1,z) = e- z/2 z 1/2+t *M{l/2 + p - K,l + 2p,z) (51) 
W(k, fi, z) = e~ z/2 z 1/2+ ^U{l/2 + p - k, 1 + 2/jl, z) (52) 

Taking into account the singularity at £ = it is required 
that the solution is bounded and well-behaved as £ — > 0. 
Following Abramowitz and Stegun (|1965p for small ar- 
guments z = 2i/3o£/°o the confluent hypergeometric func- 
tion of the second kind, U cx In z — > —00 as £ — > implying 
is an unsuitable solution. However as \z\ — > the func- 
tion .M(a,i>, 0) = 1 provided b is not a negative integer, 
in this case b = 1 implying M is a bounded, well-behaved 
solution. Hence we choose 



5(tT) = C 2 C 1/2 e 



M(- 



iE \0,*M) (53) 



Co/?0 ' 



Rewriting M as a series of modified Bessel functions of 
the first kind, I, 



M(k, p, z) 



e - Z /2 z l/2+ M x 



n=0 



D n (l/2 + p- K,l + 2p)I n (z) 



(54) 



and using the asymptotic expansion of /, the solution 
takes the form 



^,r) = (2<)- 1 / 2 C 4 exp [iufc/ 



0-0 



iPot/vo] (55) 



where D is a set of real and imaginary coefficients and 
C4 is an arbitrary constant. In the uniform equilibrium 
situation the electrostatic oscillation behaves as described 



in section 2.1 at lo 2 



2uj 2 



j 2 (Diver et al. l2"0"D2")k in 



dynamical equilibrium the electrostatic oscillation is being 
convected at the plasma flow speed oq. 

3.2. Convective solution 

If p = S = and 9 ~ 1/f ^ 9{t) this yields 



E 
A 
a 


( 
\ 
P 



-(S cr + Scr )' 

-(SoA)' 



p + /?ox 


-CTOX' - CToX/C 

-2A ( r /e 



(56) 
(57) 
(58) 
(59) 
(60) 
(61) 
(62) 



Rearranging equation (f6"2"|) for A and substituting into ([571) 
gives a differential expression for p. The same result can 
be achieved by rearranging (|59[) for x an d substituting 
into (|61|). Additionally integrating (|56|) with respect to r 
and using (|58| gives a differential equation describing the 
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evolution of E. Therefore, 



E = -<t E' + C 5 (63) 

A = -cr A' (64) 

X=-vo(Zx)'/Z (65) 

p = -M^)7e (66) 

where C5 is a constant. These have the general solution, 

£(£,t) = f(T-t/ot>)-tC B /* (67) 

A(^,r) = /(r-^/ao) (68) 

= i/(r-e/a ) (69) 



where G = p,x and / is an arbitrary functions. Here E, G 
and A are being convected at the streaming velocity of the 
plasma <jq. Assuming G, A, E ~ exp (— iwr) gives the par- 
ticular solutions oc exp [— iu>(r — £/°"o)] describing a per- 
turbation in the plasma variables propagating through the 
plasma at the streaming velocity <tq. 

3.3. General solution 

If we prescribe 6 = /(7)/£ where 7 = fc£ — (a; ± ko~o)T = 
fc£ — ilr, looking at (|43| this implies /? = 77/(7) /£, where 77 
is a constant. Substituting this expression into (|42|) yields 

= (H - pkrj) / 7 pry/ 
^ 2E 2E £ 1 ' 

Substituting ([70)) into (|3"Tj) and further prescribing a = 
A/ (7)/^, where 77 is a constant, requires 

(aofc - 0)A/C = (fi - pfc7 ? )/3 /2E (71) 
-cr A/£ = pM^o (72) 

for cr, £, and /3 to be consistent. Substituting the ex- 
pressions for cr, C, 9 and /3 into equation (|4"0)) defines the 
function / by 

f 

A (cr fc - 0) / 77 + [A(T + B (cr fc - O)] ^ 



- (1 - - foX) I = (73) 

where 

A = (0 - pkn) /2E (74) 

B = jV2E (75) 



The solution of which is a Bessel function of non-integer 
order which shows presence of Dopplcr effect as electro- 
magnetic wave propagates in the plasma that is streaming 
in the radial direction with a velocity o~q. 



4. Discussion and Further Developments 

The linear perturbation analysis presented here has shown 
some of the dynamical responses of the magnetised 
electron-positron streaming plasma. These offer guidance 
in respect of possible wave modes and plasma stability. In 
the linear regime the electrostatic, convective and general 
solutions describe the simplified streaming system per- 
turbed from equilibrium. The non-linear system of equa- 
tions does not have a closed-form analytical solution so 
the problem has to be solved numerically, work on which 
is currently still in progress. The linear analysis presented 
here will help in the development of the numerical simu- 
lations. Future considerations include extending the cold 
plasma treatment to a kinetic one. In this context the pos- 
sibility exists of coupling Bernstein modes (Laing & Diver 
I2005|) (electrostatic waves) to electromagnetic modes. 

These results imply a change in the study of pulsar 
radiation mechanisms. In pulsar magnetospheres the elec- 
tromagnetic field distribution of the star is the superposi- 
tion of the underlying dipolar electromagnetic field of the 
star, plus the self- field of the flowing plasma (da Costa 
& Kahn 111)521 da Costa et al. I200ip . Collective processes 
in the pulsar rest frame depend very strongly on the lo- 
cal plasma and field conditions. Earlier work (Diver et 
al. I2002p explored quasilincar coupling processes in a sta- 
tionary inhomogencous plasma. This article presents the 
preliminary analysis of the stability of the full non-linear 
mechanism in a streaming plasma prior to a full numerical 
simulation. 
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